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We perform a nonperturbative determination of the 0(a)-improvement coefficient csw and the 
critical hopping parameter Kc for Nf=3, 2, flavor QCD with the RG-improved gauge action using 
the Schrodinger functional method. In order to interpolate csw and Hc as a function of the bare 
coupling, a wide range of /3 from the weak coupling region to the moderately strong coupling points 
used in large-scale simulations is studied. Corrections at finite lattice size of 0{a/L) turned out to 
be large for the RG-improved gauge action, and hence we make the determination at a size fixed 
in physical units using a modified improvement condition. This enables us to avoid 0{a) scaling 



is used in numerical simulations 



' violations which would remain in physical observables if csw determined for a fixed lattice size L/a 

in ' 
o 



I. INTRODUCTION 

D . .... . n 

, Fully unquenched simulations of QCD with dynamical up, down and strange quarks have become feasible [IJ 

thanks to the recent development of algorithms and computational facilities. However, it is still very de- 

. , I _ manding to control discretization errors below a few percent level in dynamical QCD simulations. Thus highly 

I improved lattice actions are desirable to accelerate the approach to the continuum limit. 
' 

' The on-shell improvement of the Wilson quark action through 0{a) requires only a single additional term, 

i.e. the Sheikholeslami-Wohlert (SW) term In Ref. we determined csw in three-flavor QCD for the 
plaquette gauge action, using the Schrodinger functional method 0, S S 111- Applications of the resulting 
0(a) improved Wilson-clover quark action in combination with the plaquette gauge action suffer from a serious 
problem, however, since it was found in Ref. l3|] that this action combination exhibits an unphysical first-order 
phase transition at zero temperature in the strong coupling regime (/3 < 5.0). 

We also found in Ref. that such a phase transition weakens, and possibly disappears, when the gauge 
action is improved. In this work, motivated by this observation, we extend the determination of csw for the 
case of the RG-improved action [l^ for gluons for Nf=3, 2, flavor QCD. 

We explore a wide range of f3 to work out the interpolation formula as a function of the bare coupling. The 
critical hopping parameter Kc in the O(o)-improved theory is also obtained. 

In the Schrodinger functional method, csw is determined such that the axial Ward-Takahashi identity is 
satisfied for a given finite volume. Since the linear extent L of a finite lattice provides an energy scale a 
determination of csw generally involves corrections of order a/ L. We find that this correction is sizable for the 
RG improved gauge action. If the determination of csw is made for a fixed value of L/a, observables calculated 
in subsequent simulations using such csw would suffer from 0{a) scaling violations. To avoid this problem, we 
modify the standard improvement condition and determine csw at a fixed physical size L. Similar considerations 
have been made in the determinations of some other 0{a) improvement coefficients in Ref. [TTl Il2j|. 
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This paper is organized as follows. In Sec|nl we briefly recall the Schrodinger functional method, mainly 
to fix notations. In Sec. IIIII corrections at finite lattice size of 0{a/L) that affect csw are discussed, and our 
modified method and one-loop calculations relevant for the subsequent analyses are given. Section lTVl is devoted 
to describing our numerical results, and Sec. to systematic uncertainties in them. Our conclusions are given 
in Sec. lVIl A preliminary report of this work has been made in Ref. (ISj . 



II. SCHRODINGER FUNCTIONAL METHOD FOR THE DETERMINATION OF csw 



We briefly introduce the setup of the Schrodinger functional (SF) method and the improvement condition 
developed in Refs. HIEIlli- 



A. SF setup 



Consider the SF defined on a four dimensional hypercubic lattice with a volume L'^ x T and the cylindrical 
geometry, i.e., the periodic boundary condition is imposed in the spatial directions and the Dirichlet one in the 
temporal direction for both gauge and quark fields. At the temporal boundaries xq = and T, the following 
conditions are imposed on the link variables and the quark fields: the spatial link variables on the boundaries 
are fixed to the diagonal, constant SU{3) matrices given by 

t/fc(x,xo)|^^ = exp[aCfc] , ;7fe(x, xo)!^^^^ = exp [aC^] , (1) 

. / -5 \ 
> - 7^ I 2 I , (2) 

while all quark fields on the boundaries are set to zero. 
We use the RG-improved gauge action given by, 

2 
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^5 = 4 X [ E <-(^o) (1 - P^^A^)) + E <(^o) Tr (l - Ri]i^Hx)) ] , (3) 



f 1 X 2) 

where P^,,y(a;) denotes a 1x1 Wilson loop on the ji-u plane starting and ending at x, and Rli.u {x) a 1x2 
rectangular loop with the side of length 2 in the v direction. These terms are added up with proper weights, 
w^j^(xo) and wj}jj{xo), respectively. In ordinary simulations with the periodic boundary condition in the temporal 
direction, the weights are given by wj^j^—3.6A8 and w^^^— 0.331 independently of xq. In the SF, these weights 
are modified. Among several possible choices, we select the choice B defined in Ref. 14] in this work, 

, - X (3.648) at t = or T and /i, ^4 
w^A^o) = S 2 , (4) 

3.648 otherwise 

at t = or T and fi, v ^4 

w^ui^o) = ^ X (-0.331) at i = or T and ^=4 • (5) 
—0.331 otherwise 

The 0(a)-improved Wilson quark action ^ is given by 

Sq = ^qxD^yQy, (6) 

x,y 
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with the field strength tensor Fx.fj,i, defined by 

F..^. = I {{P^,..{x) + P.,-,.{x) + P-^,-.{x) + P-,,f,{x)) - (h.c.)} , (8) 

o 

and a^,^ = {i/2) [7^,7,/]. The last term in Eq. (0 is the only counter term to get rid of 0(a) errors present 
for on-shell quantities on the lattice. At tree level, csw=l- For the 0(a)-improvement of the SF, we need to 
add extra terms made of the gauge and quark fields at boundaries to the lattice action. However, since these 
counter terms affect the PCAC relation used in the following calculations only at 0{a^) or higher, they are not 
necessary for the determination of cgw- 



B. PCAC relation 



We determine csw by imposing the PCAC relation 

i {d^ + d;) ^r.„p,^ = 2m,P^ (9) 

up to 0{a^) corrections. The pseudo-scalar density operator, axial vector current and its 0(a)-improved version 
are given by 

P- = ^Pj^r'i^, (10) 

= i'l^l^r'^i^, (11) 

^Lp,^ = + + (12) 

where 9^ and d* are the forward and backward lattice derivatives, and denotes the generator of SU{Nf) 
flavor symmetry acting on the flavor indices of the quark fields 'ip and ip- 



We measure two correlation functions, 



fA{xo) = -^^^(AS(x)O"), (13) 



'/ 

fpixo) = -j^j-j{P''{x)0'^), (14) 

where x = (a;o,x), and (• • ■ ) represents the expectation value after taking trace over color and spinor indices 
and summing over spatial coordinate x. The source operator is given by 

= a6^C(y)75r"C(z), (15) 

C(X) = C(X) = J^y (16) 

(5p(x) Sp{x) 

where p(x) is the quark field at a;o = and is set to zero in the calculation of Ja and fp. The bare PCAC quark 
mass is then calculated using and fp through the PCAC relation Eq. © as 

m{xo) = r{xo) + CAs{xo) (17) 
r{xo) - \{do + d*o)fA{xo)/fp{xo) (18) 

sixo) = ^adodQfp{xQ)/fp{xo). (19) 

Using the source operator on the other boundary 

^ a«^C'(y)75r'^C'(z), (20) 
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where (' is the boundary field at xq = T, we can calculate another set of quantities m'(a;o), r'{xo) and s'{xd) 
from the correlation functions defined by 

UT-xo) = (21) 
f'p{T-xo) = -j;^{P-{x)0''''), (22) 

A naive improvement condition would be m(a;o) =iti'{xo). However, this condition requires a nonperturbative 
tuning of ca as well as of csw- To eliminate ca from the determination, it was proposed in Ref. Q to use an 
alternative definition of the quark mass given by 

M(xo,yo) = m[xo) — — — s(xo), (23) 

s[yo) - s'[yo) 

M (xQ^ya) = m xo -— --—s{xo). 24 

s (yo) - s{yo) 

with which csw is obtained at the point where the mass difference 

AM(xo, yo) = M {xo, ya) - M'{xo, yo) (25) 

vanishes. In principle, we can take an arbitrary choice for (xo,?/o)i since different choices result only in O(a^) 
differences in physical observables. We follow the ALPHA Collaboration and use (a;o,yo) — (3T/4,T/4) for 
AM, and (r/2,r/4) for M. In the following, M and AM without arguments denote M(T/2,T/4) and 
AM(3T/4,T/4), respectively. 

In previous studies, csw has been determined through the conditions 

M{glL/a) = 0, , . 

AM{glL/a) = AM(0,L/a), ^^^^ 

at a given and L/a. AM(0, L/a) on the right hand side, which is the tree-level value of AM((7q, L/a) at the 
massless point, is necessary in order that the resulting csw reproduces its tree-level value (csw=l) in the weak 
coupling limit. In the next section, we address the issue of corrections at finite lattice size, and propose a new 
condition to avoid the problem. 



III. CORRECTIONS AT FINITE LATTICE SIZE AND MODIFIED IMPROVEMENT 

CONDITIONS 



A. corrections at finite lattice size 



In the standard approach, we first calculate M{gQ, L/a) and AM{gQ, L/a) for a set of values of csw and n. 
The results are fitted as a function of csw and k to find csw(ffo' ^1^) ^nd Kc{gQ,L/a) satisfying Eq. at a 
given value of g^ and L/a. The asymptotic a dependence of cswi9o,L/a) and Kc(gQ, L/a) obtained in such a 
way is expected to be 

csw{9o,L/a) = csw(.9o:Oo) +CL • (a/i) + CA • (aAqcD) +0((a/i)^(a^AQCD/i),(aAQCD)^), (27) 
Kc{gl,L/a) = Kc{gl,oo) + kL-{a/L) + kA-{aAQCD)+0{{a/L)\{a^AQCD/L),{aAQCD)^), (28) 

where cl, ca, and k\ are unknown coefficients. (In practice, a logarithmic dependence on a/ L also appears, 
but it does not alter the following discussion, and hence not written explicitly.) 

Consider an on-shell physical quantity Q, and let Q'*"(a) be the value obtained on a lattice with lattice spacing 
a using the SW quark action with a choice of the improvement coefficient Cg-y^. We expect the discrepancy 
between Q and (5''**'(a) in the measured value to be 

Q - Q^-*\a) = q • (c|™ - cswiol^)) ■ (oAqcd) + 0(a^ A|cd), (29) 



5 



where q is an unknown constant assumed to be 0{\). Hence, if one uses Cg{y = csw{9o,'X)) in the simulation, 
the 0{a) error is absent, while if one uses csw(5oi L/a) in Eq. (|27() . the above expression results in 

Q - Q'-^\a) ^q-CL- (a/L) ■ (aAqcD) + 0{a^ A^cd) + 0[aKciMa/Lf). (30) 

While the scaling violation appears to start from 0{a?), it is actually linear in the lattice spacing if one deter- 
mines csw(5o' L/'^) with a fixed value of a/ L. Indeed, previous studies determining csw have used certain fixed 
values of -L/a, e.g. 8, independently of 

In Ref. H, we studied the magnitude of the corrections at finite lattice size in csw for the plaquette gauge 
action. The coefficient defined in Eq. (|27|l was evaluated in one-loop perturbation theory in the same SF 
setup, and it was found that the effect on csw does not exceed 3% when L/a— 8 for (3 > 5.2. We have repeated 
the same perturbative analysis with the RG-improved action, and observed a sizable effect of about 15% at 
/3=1.9, around which large-scale simulations are carried out. This enhancement of the one-loop correction for 
the RG improved action is mainly due to the larger value of the bare coupling compared to that for the plaquette 
gauge action for realizing the same value of the lattice spacing. 



B. modified improvement condition 



2 r - n (31) 



We propose to resolve the problem due to the sizable corrections explained above by introducing a fixed 
physical length L*, and determining csw at the fixed physical volume L*^ x T* (T* = 2L*). If one uses csw 
thus determined, L in (|30|) is replaced by L* and scaling violations are O(a^). 

The actual procedure we use runs as follows. Instead of Eq. H26I) . we impose a modified improvement condition 
given by 

MiglL/a) = 
AM{glL/a) = 

to calculate csw(5oi L/a) and Kc{9o; L/a). The results are converted to csw(5o' L* / a) and Kc(3o, L* / a). To do 
so, we must know the value of L* / a or l/a at that value of ^q, which we obtain through the two- loop (3 function, 

aAi = exp (-^) {b,gl)-'^/^'\ (32) 
47r 2 \ Z " 3 



^0 = TTTW ( - o^/ ) ' (33) 



^1 = A-A-NI~nA'^N^^1-\\. (34) 
(4^)4 \ ?, ^ ^ V 3 Nj) ^ ' 

The transformation from csw(5o' L/a) and Kc(5oi L/a) to those at L* / a are made through 

csw{9f),L* /a) = csw{9l,L/a) + 6csw{go^L/a;L*/a), (35) 

i^c{9l,L*/a) = Kc{go,L/a) + dKc{gl,L/a;L*/a), (36) 

where 

Scsw{9o,L/a;L*/a) = -cl^ig^, L/a) + c^^igl, L* /a), (37) 

Sn,{glL/a;L*/a) = -t^J {gl L / a) + t^J {gl L* / a), (38) 

and Cg^((7o, L/a) and t^'^lg^, L/a) are calculated at the one- loop level for the same SF setup at the given value 
of L/a. 

It turned out that the tree and the one-loop coefficients for csw and Kc have a significant a/L dependence. 
To describe this dependence precisely we fit them to a Fade or a polynomial-like function of a/L as 

SO) (T i„\ - 1 + fli [a/L) + aa ja/L)'^ + 03 (a/L)^ 

1 + 61 [a/L) 

4w(^/«) = 0.113 +(ci-diln(L/a)) (a/L) + (c2-d2ln(L/a))(a/L)^ (40) 

4*') (L/a) = i + fci (a/L) + A:2 (a/L)2 + fcg (a/L)3 + fc4 (a/L)^ (41) 
8 

(L/a) = 0.002760894+ (^1 -mi ln(L/a)) (a/L) (^2 -m2ln(L/a)) (a/L)2. (42) 



(0)/;-, X _ -L ^ ui yu./ -tv; ^ c^2 V"./ -t^; -ri^av^/^y /oQ\ 
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The coefficients are given in Tabled We note that the one-loop coefficients have an Nf dependence due to the 
tadpole diagram, although it vanishes in the large volume limit. 

In our actual determination, we define L* by L*/a = 6 at /3 = 1.9, L*/a = 6 at /? = 2.0 and L*/a = 6 at 
P = 2.6 for Nf=3, 2, flavor QCD, respectively. In Table HTHlVI numerical values oi P — Q/qq, L/a and L* /a in 
our simulations for Nf=3, 2, and cases are summarized. In these tables, we also show the numerical values of 
5csw{9o,L/a]L* /a) and SKcig^, L/a; L* /a). For large values of /3, the perturbative corrections are small and 
hence reliable. On the other hand, if L/a are close to L* /a, the corrections needed for the conversion from L 
to L* should again be small. Since we fix L* at strong coupling, the corrections, Eqs. H37I) and l|38|) . are small 
at both ends of our range of (3 as one can see in the Tables. 



IV. NUMERICAL SIMULATIONS 



A. parameters and algorithm 

The numerical simulations are performed with Nf—3, 2 and degenerate dynamical quarks on a [L/a^ x 
2{L/a) {L/a— 8 or 6) lattice for a wide range of (3. The simulation parameters are summarized in Tabs. HTIlIVI 
for Nf—3, 2, 0, respectively. 

We employ the symmetric even-odd preconditioning introduced in Refs. [isl for the quark matrix D. 
Calculation of D^^ is made with the BiCGStab algorithm with the tolerance parameter ||i?i||/||i3|| < 10~^^, 
where Ri — DXi — B is the residual vector and Xi is an estimate for the solution X in the i-th BiCGStab 
iteration. 

We adopt the standard HMC algorith m IITII for the Nf—2 and flavor cases. For the three-flavor case, the 
polynomial HMC (PHMC) algorithm [TgL I18| is applied to describe the third flavor, employing the Chebyshev 
polynomial P[D] to approximate D~^. In order to make the PHMC algorithm exact, the correction factor 
-PcoiT = det[VF[£']] with I^[D] = P[D]£) is taken into account by the noisy Metropolis method 0. The square 
root of ^^[Z?], which is required in the Metropolis test, is evaluated with an accuracy of 10""'^^ using a Taylor 
expansion The order of the polynomial A'poiy is chosen so that an acceptance rate of about 70% 

or higher is achieved for the Metropolis test. 

In the calculations of aM and aAM, fx and {X=A or P) are first evaluated at every trajectory, and they 
are combined to produce aM and aAM. The bin size dependence of the jackknife error of aM is investigated 
in the range iVbin = l--/Vtraj/20. We adopt A^bin giving the maximum error in this range in the error analyses 
in the following. 



B. results 

The trial values of csw and k at which simulations are made are summarized in Tables. IVHVIII for Nf—3, 2, 
and 0, respectively, together with the results for aM and aAM and the number of trajectories accumulated. 
In order to obtain csw(5o' ^f'^) '^c(ffo' ^/^) satisfying Eq. at each /3, we make fits of those data using 
the functional forms, 

,(1) ^(2) , 

aM = a^, + + JM_ + 4Vcsw + 4/ csw' + — csw, (43) 

aAM = a^M + ^ + ^ + ^sw + 41 4^ + ^ csw (44) 

The results for csw{9Q,L/a) and Kc{go,L/a) obtained with the fits, and the adopted functional form are 
tabulated in Tabs. IVIIIhlXl The details of the fit procedure are as follows. In Figs. ^3 we plot data on the 
(aM, aAM) plane for Nf—3, 2,0, respectively. For those data for which the origin (0, 0) is contained in or close 
to the data region, we make a fit leaving only the constant and linear terms in Eqs. (|43|l and H44|l . This applies 
to all cases except for the three- fiavor simulations at (3 < 2.2, and the dotted lines in the figures show the fit 
results. 



7 



In the three-flavor simulations at /? < 2.2, the region of negative aM is not covered, and the origin is missed 
by the data. This happens because the PHMC algorithm tends to fail at vanishing or negative PCAC quark 
masses at low /3 due to large quantum fluctuations. Thus, at /? < 2.2, we are forced to extrapolate the data. 
In the extrapolation, three functional forms are examined: (i) linear, (ii) quadratic without the cross terms, 
and (iii) quadratic with the cross terms. At /3=2.20 and 2.10, a linear function well fits the data, and we take 
this in the following analysis. The data at j3—2.QQ and 1.90 require the quadratic term, but it turns out that 
including the cross terms does not reduce x^/dof significantly from that without the cross terms, and leads to 
csw(5oi -^1 ^) '^c(.9oi ^/^) consistent within one standard deviation. Thus, we adopt the quadratic function 
without the cross terms at these /3, and (Im and d/s^M are always set to zero throughout this analysis. 

Next, csw{gQ,L/a) and Kc(5§,i/a) are transformed into those for the desired lattice volume, {L*/a)^ x 
2{L*/a), along the line presented in Sec. 1111 lH Using Eqs. H35|) . H3t)|) and the 6csw and 6kc given in Tables [TTI 
HVl we obtain cswigl.L* /a) and Kc[gl,L* /a) shown in Tables IXll - IXIllI Notice that in Table IXll there are 
three results for /3=2.0. The first and second one are obtained by transforming the data with 8^ x 16 and 
6^ X 12 to those for L* /a ~6.805, respectively, and the third one is obtained by simply interpolating the two 
raw values at L/a = 8 and 6 in Table IVnTI to L* /a ~6.805, for which the corrections at finite lattice size are 
essentially corrected nonperturbatively. The two raw values, 1.670(56) at L/a — 8 and 1.632(45) at L/a ~ 6, 
are very close to each other and consistent within the error, and hence the linear interpolation to L* /a 6.805 
is more reliable than the perturbative procedure. Similar observations are made at the second smallest P in 
each Nf flavor simulation, namely at /3=2.10 for Nf = 2 and at /?=2.70 for Nf = 0. Thus, at these (3 the 
result interpolated to L* /a is adopted as our final result, and used in the following analysis. At the same time, 
it is worth noting that in all three cases the one-loop corrections have the right sign, which indicates that the 
one-loop correction dominates over higher loop corrections. Furthermore, the discrepancy between the results 
corrected perturbatively and nonperturbatively is found to be 5%, 3% and less than 1% for the Nf=3, 2 and 
cases, respectively, while the size of one-loop correction itself at these (3 is 6-7%, 5% and 2-3%. From this 
observation, we expect that the size of the one-loop correction gives a conservative estimate for the unknown 
higher loop corrections for all f3. 



C. interpolation formula 

Our final results for csw(5o' ^*) ^ function of gg are shown in Fig.^jfor Nf = 3, 2, fiavor QCD. When we 
interpolate csw, not all available data are used in the fit. As mentioned in Sec. IIII 51 the corrections at finite 
lattice size estimated perturbatively is small only around the high and low ends of /3 due to our choice of L*, 
while in the middle range corrections may be significant. Therefore, we use data only if the correction is less 
than 5%. In the three flavor case, the data at /3=12.0, 8.85, 2.2, 2.1, 2.0, 1.9 are employed. As a consequence, 
we obtain the foUowings interpolation formula, 

csw(go, = 1 + 0.113 gl + 0.0209(72) {g^f + 0.0047(27) (.9^)3, {^^/doi = 0.58). (45) 

For Kc shown in Fig. [S] the corrections are smaller than 5% for all value of (3. Including all data in the fit we 
obtain 

Kc{gl,L*) = 1/8 + 0.003681192 g2 +0.000211(43) (g^)2 

-1-0.000067(66) (.g^)^ - 0.000038(21) (gl)^. {x^/cM = 1.1) (46) 

When performing the above fits, the tree and one-loop coefficients are fixed to the perturbative values at infinite 
volume. This is justified since, as seen in Table UTI L* /a grows very rapidly with /3, and hence a/L* corrections 
in Eqs. H39|) - (|42|l are all negligibly small near the continuum limit. We also note that the tree and one-loop 
coefficients in the infinite volume limit do not depend on iV^ , and hence the same values are used in the analysis 
for the Nf = 2 and cases given below. 

The interpolation formula for csw in two-flavor QCD is calculated in the same fashion as in the three-flavor 
case. In this case, the sizes of the correction at finite lattice size are acceptable (< 5%) at /3=12.0, 8.85, 5,0, 
2.2, 2.1, 2.0. We first try a polynomial form as before, and obtain 

cswiglL*) = 1 + 0.113 <?2 + 0.0158(63) (5o')' + 0.0088(24) (g2)3, (xVdof = 4.68), (47) 

which is denoted by a dashed line in Fig. ^ A sharp rise of the data points near ^q^S.O is not described well 
by this polynomial form, while in the three-fiavor case the polynomial worked well over the whole range of (3 
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we studied. An alternative is a Pade function, with which we obtain 

, ,„2 _ 1 - 0.212(9) ffg - 0.0108(38) {glf ~ 0.0083(19) {gjf 

cswi9o,L )^ 1 - 0.325(9) g§ ' (x/doi = 2.n), (48) 

This fit, denoted by solid line in the middle panel of Fig. 0] interpolates our data very well. Since this formula 
has a pole at §^=3.08(8), its use is restricted to (7o^3.0. For Kc, we use all available data to obtain 

Kc{go,L*) = 1/8 + 0.003681192 5^ + 0.000227(58) (5^)2 

+0.000093(84) (.9^)3 - 0.000049(24) (.9^)^ (xVdof = 0.98), (49) 

for a polynomial, and 

2 1/8 - 0.0356(23) g§ - 0.00089(8) (g§)^ - 0.00009(6) (gg)=^ 2/, 

M5o,^)- 1- 0.314(18) g2 ' (X /dot - 0.35), (50) 

for a Pade function. These results appear in the middle panel of Fig. as dashed and solid line, respectively. 
It is interesting that the pole positions for csw and Kc are consistent with each other. This seems to indicate 
that above gp ~3.0 the Wilson quark action cannot be improved in this fashion consistently for the Nf=2 case. 
All in all the Pade fits provide a more satisfactory interpolation of the Nf=2 data, and we take them as the 
main result for the Nf=2 case. We have also applied a Pade function for csw in the Nf=3 case. However, in 
this case the resulting fit lies on top of that for a polynomial over the range of (3 we used, and the position of 
pole can be determined only poorly. Hence there seems no reason to favor the Pade fit over the polynomial for 
interpolating the data. The difference between the Nf—2 and 3 cases probably arise from the fact, empirically 
known, that the Nf—2 lattice is coarser than the Nf=3 lattice at the same value of gp- Indeed, a sharp rise of 
improvement coefficients was previously seen for the plaquette gauge action toward coarse lattices 0, 1^ • 

In quenched QCD, the size of the correction is smaller than 5% for all availble data, and we use all data to 
obtain 

csw(5o:^*) = 1 + 0.113^2 + 0.0371(54) (g2)2_ 0.0036(26) (5^)3, (^Vdof = 4.09), (51) 
Kc{gl,L*) = 1/8 + 0.003681192^^ + 0.000293(37) (gg)2 

-0.000053(65) {g^f + 0.000008(24) (5g)^ (xVdof = 0.46). (52) 

In Rcf. '2^, the authors performed a one-loop determination of Cg"!^ with conventional perturbation theory, 

and reported a very precise value Cg\Y=0. 11300591(1) in the infinite volume limit. Changes in our results due 
to the use of this value in above analyses are expected to be negligibly small. 



V. SYSTEMATIC ERRORS 



There are two sources of systematic errors in our analysis, both related to the conversion to a fixed physical 
length scale L*, one being the use of the two- loop (3 function to estimate L* as a function of gQ, and the second 
being the use of one-loop perturbation theory for correcting the value of csw from L to L*. 

In order to examine the magnitude of uncertainties from the first error, we go through the analysis using 
the three-loop f3 function. Since the three-loop term of the lattice /? function is not available for the RG- 
improved gauge action, we take the value for the plaquette gauge action. Thus the following argument is only 
semi-quantitatively valid. In this case, Eq. 1)32(1 is replaced with 

uAl = exp (-^) (bogl)-'"/^''" X (1 + qgl) , (53) 

where g=0.18960350(l), 0.4529(1), and 0.6138(2) for A^/=0, 2, and 3 [2l |, resp ectively. With this function, we 
estimate L*/a, Scsw, and Skc with Nf=3, which are tabulated in Table IXIVl Comparing with Table ITTI it is 
found that L* /a changes significantly while the changes in dcsw and Suc are at most a few percent and hence 
small. Thus we conclude that the uncertainty from scaling violation in the lattice spacing is negligible. 



9 



In order to discuss the uncertainty of one-loop corrections, we write csw(5oi^*) determined through our 
procedure as 

+gl (c(^)(a/L) - c^'Ha/L*)) + 0{g',). (54) 

In other words, Eq. H54|l represents the difference between csw(5o' ^*) csw(5o' terms of perturbative 

series with coefficients c''*^(a/L), where c^-^^{a/L) vanishes as L ^ oo. Since we have corrected the mismatch 
between csw{gQ,L/a) and cswig^iL* /a) only at the tree- and one-loop level, the unwanted a/L dependence 
remains at two-loop and higher. Replacing Cg^ in Eq. (|29|l with Eq. H54|l . we obtain 

Q - Qi-"(a) = (^c^^\a/L*)+g'oC^'Ha/L*)+gtcS'Ha/L*)^ • (aAqcd) 

+.g4 (c(2)(a/L) - c(2)(a/L*)) • (oAqcd) + ©(so'aAqcD a/L) + 0{a^ A^cd), (55) 

where we omit an unknown 0(1) overall coefficient g, because it is not relevant in the following discussion. If 
you expand c^-^\a/L'-*^) around a/L'--*^ — 0, the first term in Eq. H55|l behaves ^ o^Aqcd/L* ^ O(a^) because 
L* is fixed. The second term behaves like ~ 5o(a/L — a/L*)(aAQCD), which gives 0{a) scaling violation because 
a/L is fixed. As a results, the leading scaHng violation could be 0{a) rather than O(a^). However it should be 
emphasized that when we obtain the interpolation formula we only used the weak coupling and strong coupling 
regions because in these regions the perturbative errors are expected to be under control for the following 
reasons. In the weak coupling region, L/a and L* /a are different by several orders of magnitude, but the 
couphng is very small, and hence the size of 0{g'^{a/L — a/L*)(aAQCD)) is expected to be as small as the size 
of the one- loop corrections. On the other hand, in the strong coupling region, L/a and L* /a are close to each 
other, and again the remaining 0{a) scaling violation, 0{g*{a/L — a/L*)(aAQCD)), should be small. We also 
saw in Sec. lIVBl that the size of perturbative errors is roughly the same as that of the one- loop correction itself. 

Most importantly, at our strongest and the second strongest couplings around which large-scale simulations 
are performed, there are no perturbative errors in csw due to our choice of L* and interpolation to L* at the 
second strongest couplings. Thus we believe 0(a) scaling violations are well below O(a^), though we need to 
check this in future work. 



VI. CONCLUSION 



In this work, we have performed a nonperturbative determination of the 0(a)-improvement coefficient cgw of 
the Wilson quark action with the RG-improved gauge action for Nf—3, 2, and flavor QCD. The corrections at 
finite lattice size turn out to be sizable, and are taken into account by modifying the improvement condition and 
carrying out the determination at a fixed physical length scale of L*. While we have to resort to perturbation 
theory to incorporate the corrections, we have attempted to choose L* at a moderately strong coupling, close 
to the range of lattice sizes of order ~ 2 GeV where physics simulations are practically made, so that their 
magnitude are reasonably under control. 

Using the data for csw thus obtained over a wide range of we have determined the interpolation formulas, 
given in Eqs. (|45|l . H48I) and (|51|l . which represent the main results of this work. These results do depend on L* 
chosen, but the removal of 0{a) scaling violations in physical observables hold independent of the value of L*. 

As a byproduct, we have also obtained the interpolation formula for Kc, Eqs. H46|l . H5U|I and H52I) . which may 
be useful to locate simulation points. 

The three-flavor results reported here are already being used in a large-scale simulation aiming to carry out 
a systematic evaluation of hadronic observables for the realistic quark spectrum incorporating the dynamical 
up, down and strange quarks. The preliminary results have been reported in Refs. p^ . 
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TABLE I: Finite-size coefficients in Eqs. lO^-lH^. 



'-sw 




^/=0 




^/=3 


ai 


-3.4415 


Cl 


-4.5736 


-6.2641 


-7.1094 


ai 


-5.0248 


C2 


-3.3402 


-8.0488 


-10.403 


as 


11.1475 


rfi 


-1.1681 


-1.5466 


-1.7359 


61 


-3.9702 


rf2 


-8.9448 


-14.306 


-16.987 














iV/=0 


Nf = 2 






0.260982x10"'' 


/i 


0.101302x10"^ 


-0.224650x10"^ 


-0.387626x10"^ 


fc2 


-0.845333x10"^ 


/2 


0.162496x10"^ 


0.862878x10"^ 


0.481835x10"^ 


fc3 


-0.103610x10"^ 


mi 


0.547826x10"^ 


-0.507665x10"^ 


-0.155835x10"=* 


^4 


0.751742x10"^ 


m2 


0.882220x10"^ 


-0.136413x10"^ 


-0.645729x10"^ 



TABLE II: Inverse coupling [3 and lattice size L/a chosen for the three-flavor QCD simulation. L* /a is estimated by the 
two-loop f3 function assuming L* /a = 6 at p—1.9. Finite-size corrections Scsw and 5kc calculated with Eqs. (1371 and 13811 
are also shown. 





L/a 


L'/a 


<5csw(5'o, L/a; L* /a) 


Sucigl, L/a; L* /a) 


12.00 


8 


7.51x10" 


5.51x10"^ 


6.35x10"^ 


8.85 


8 


8.46x10" 


1.42x10"^ 


7.95x10"^^ 


5.00 


8 


3.81x10^ 


5.14x10"^ 


1.23x10"* 


3.00 


8 


2.50x10^ 


1.14x10"' 


6.80x10"^ 


2.60 


8 


1.48x10^ 


1.08x10"' 


1.34x10"^^ 


2.40 


8 


1.14x10^ 


8.70x10"^ 


-8.82x10"'^ 


2.20 


8 


8.78 


3.42x10"^ 


-9.84x10"'^ 


2.10 


8 


7.73 


-1.59x10"^ 


5.70x10"^ 


2.00 


8 


6.81 


-9.36x10"^ 


3.85x10"^^ 


2.00 


6 


6.81 


1.10x10"' 


-5.08x10"=^ 


1.90 


6 


6 












TABLE III: Same as 


Table HH but for two-flavor QCD. 




/3 


L/a 


L'/a 


Scsw{L/a; L* /a) 


5ti{L/a; L* /a) 


12.00 


8 


2.35 xlO'^ 


2.43x10"^ 


5.93x10"^ 


8.85 


8 


3.66x10'' 


1.01x10"^ 


7.38x10"^ 


5.00 


8 


2.45x10^ 


4.63x10"^ 


1.08x10"" 


3.00 


8 


1.98x10^ 


9.51x10"^ 


2.69x10"^ 


2.60 


8 


1.22x10' 


7.84x10"^ 


-1.81x10"=^ 


2.20 


8 


7.58 


-2.11x10"^ 


1.41x10"^ 


2.10 


8 


6.74 


-8.24x10"^ 


6.14x10"^ 


2.10 


6 


6.74 


8.37x10"^ 


-7.04x10"^ 


2.00 


6 


6 
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TABLE IV: Same as Table [TTl but for quenched QCD. 

13 L/a L*/a 5csw{L/a; L* /a) 5k,{L / a; L* / a) 

24.00 8 3.09x10" -1.11x10"^ 3.47x10"^ 

12.00 8 2.41x10^ -3.70x10"^ 5.08x10"^ 

8.85 8 6.33x10^ 2.24x10"^ 6.18x10"^ 

5.00 8 8.04x10^ 3.80x10"^ 6.29x10"^ 

3.00 8 9.12 2.07x10"^ -2.68x10"^ 

2.70 8 6.66 -4.81x10"^ 8.79x10"^ 

2.70 6 6.66 3.98x10"^ -8.40x10"'' 

2.60 6 6 
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TABLE V: Results for aM and a AM for three-flavor QCD. The accep- 
tance rates for the MD and the noisy Metropolis test are shown together 
with the number of MD steps per trajectory and the order of the poly- 
nomial Afpoiy used in the noisy Metropolis test. The final column gives 
the number of trajetories accumulated. 



CSW K aM a AM PacdA^MpI fcorr|-/Vpolv| NtTIj 



13 = 12.00, L/a = 8 



1.00 


0.12659 


0.01235(13) 


0.00101(13) 




0.73(2) 


100 


0.983(5) 


100 


1600 




0.12676 


0.006906(91) 


0.00072(13) 




0.75 




100 


0.979 




100 


1600 




0.12693 


0.00149(13) 


0.00087(14) 




0.73 




100 


0.973 


^! 


100 


1600 




0.12709 


-0.00368(13) 


0.00092(16) 




0.75 


1) 


100 


0.970 


5) 


100 


1600 


1.05 


0.12659 


0.008565(98) 


-0.00009 


li 




0.72 


2) 


100 


0.984 


4) 


100 


1600 




0.12676 


0.00283(11) 


-0.00003 


13^ 


0.74 


1) 


100 


0.968 


6) 


100 


1600 




0.12693 


-0.00221(11) 


0.00004(17) 




0.72 


'3) 


100 


0.969 


6) 


100 


1600 




0.12709 


-0.007708(89) 


0.00023(10^ 




0.74 


2) 


100 


0.955 


5) 


100 


1600 


1.10 


0.12659 


0.00460(12) 


-0.00076 


i; 




0.72 


2) 


100 


0.981 


4) 


100 


1600 




0.12676 


-0.00097(15) 


-0.00069 


22) 


0.71 


3) 


100 


0.972 


4) 


100 


1600 




0.12693 


-0.00625(21) 


-0.00073 


16) 


0.74 




100 


0.953(9) 


100 


1600 




0.12709 


-0.01161(17) 


-0.00070 


13^ 


0.73 


■11 


100 


0.94(2) [1001 


1600 



P = 8.85, L/a = 8 



1.0141 


0.12698 


0.02316(13) 




0.00094(15) 


0.68(2) 


80 


0.989(3) 


100 


2000 




0.12730 


0.01311(11) 




0.00104 


18) 


0.66 


2) 


80 


0.988 


4) 


100 


2000 




0.12762 


0.00309(12^ 




0.00102 


15) 


0.71 


2) 


80 


0.969 


5) 


100 


2000 




0.12826 


-0.01734(10) 


0.00092 


16) 


0.70 


1) 


80 


0.943 


6) 


110 


2000 


1.0350 


0.12698 


0.02121(11) 




0.00060 


12) 


0.70 


2) 


80 


0.990 


2) 


100 


2000 




0.12730 


0.01119(16 




0.00057 


17) 


0.71 


1) 


80 


0.985 


3) 


100 


2000 




0.12762 


0.00080(13^ 




0.00060 


16) 


0.70 


2) 


'80 


0.972 


4) 


100 


2000 




0.12826 


-0.01948(13) 


0.00055 




0.71 




80 


0.937 




110 


2000 


1.0559 


0.12698 


0.01903(14) 


0.00013 


51! 


0.69 


1! 


'80 


0.987 


% 


100 


2000 




0.12730 


0.00896(11' 




0.00033 


10) 


0.68 


1) 


80 


0.975 


4) 


100 


2000 




0.12762 


-0.00130(11) 


0.00024 


20) 


0.69 


2) 


80 


0.964 


5) 


100 


2000 




0.12826 


-0.02166(13) 


0.00039 


10) 


0.70 


2) 


80 


0.928 


7) 


110 


2000 


1.0800 


0.12719 


0.00983(15) 


-0.00013(15) 


0.69 




80 


0.990 


^! 


120 


2000 




0.12753 


-0.00078(11) 


0.00007(12) 


0.68 




80 


0.989 




130 


2000 


1.1000 


0.12713 


0.00982(35) 


-0.00041 


(19) 


0.68 




80 


0.990 




110 


2000 




0.12747 


-0.00121(11) 


-0.00039(14) 


0.69 




80 


0.986 


^! 


120 


2000 



P = 5.00, L/a = 8 



1.08 


0.12958 


0.01031(33) 


0.00073(22) 


0.72 


1) 


64 


0.982 


3) 


100 


2200 




0.12974 


0.00553(17) 


0.00082(29) 


0.77 


2) 


64 


0.968 


4) 


100 


2200 




0.12989 


0.00049(35) 


0.00060(25) 


0.74 


2) 


64 


0.970 


4) 


100 


2200 




0.13004 


-0.00377(20) 


0.00070(19) 


0.74 


2) 


64 


0.962 


5) 


100 


2200 


1.13 


0.12932 


0.01027(26) 


0.00039(22) 


0.73 


1) 


'64 


0.976 


6) 


100 


2200 




0.12948 


0.00541(24) 


0.00043(27) 


0.73 


1) 


'64 


0.975 


4) 


100 


2200 




0.12963 


0.00030(24) 


0.00024(21) 


0.77 


1) 


64 


0.964 


7) 


100 


2200 




0.12978 


-0.00428(29) 


0.00042(37) 


0.74 


1) 


64 


0.950 


6) 


100 


2200 


1.18 


0.12907 


0.01002(22) 


-0.00071(20) 


0.76 


2) 


64 


0.982 


3) 


100 


2200 




0.12922 


0.00489(25) 


-0.00105(22) 


0.73 


2) 


64 


0.974 


4) 


100 


2200 




0.12937 


0.00043(43) 


-0.00074(21 


0.75 


1) 


64 


0.970 


4) 


100 


2200 




0.12952 


-0.00441(28) 


-0.00089(21^ 


0.76 


*1) 


64 


0.958 


*6) 


100 


2200 



(5 = 3.00, L/a = 8 



1.20 


0.13281 


0.02798(33) 


0.00099(26) 


0.770(9) [501 


0.988 




100 


4300 




0.13311 


0.01769(47) 


0.00036(57) 


0.77(1 




50 




0.978 


1 


100 


4000 




0.13341 


0.00813(46) 


0.00044(33) 


0.78 


1 




50 




0.960 


4) 


100 


4000 




0.13370 


-0.00036(41) 


0.00104(40) 


0.77 


1 




50 




0.937 


5) 


100 


3800 


1.25 


0.13235 


0.02774(42) 


-0.00015(61) 0.78 


1 




50 




0.986 


4) 


100 


4200 




0.13265 


0.01820(52) 


0.00003(40) 


0.76 


1 




50 




0.980 


3) 


100 


3800 




0.13294 


0.00940(38) 


0.00004(36) 


0.77 


1 




50 




0.960 


3) 


100 


4200 




0.13324 


-0.00027(40) 


-0.00029(34) 


0.773(9) [501 


0.948 


5) 


100 


3900 


1.30 


0.13190 


0.02742(65) 


-0.00052(28) 


0.76 


1 




50 




0.990 


2) 


100 


4200 




0.13219 


0.01713(64) 


-0.00040(36 


0.77 


2 




50 




0.980 


2) 


100 


4200 




0.13248 


0.00915(54) 


-0.00066(67, 


0.78 


2 




50 




0.962 


3) 


100 


3900 




0.13278 


-0.00008(41) 


0.00031(64) 


0.77 


1 




50 




0.945 


6) 


100 


4000 


1.35 


0.13145 


0.02697(35) 


-0.00098(46) 


0.77(1 




50 




0.988 




100 


4300 




0.13174 


0.01643(70) 


-0.00075(44) 


0.770(8) [501 


0.979 


5! 


100 


4100 




0.13203 


0.00800(43) 


-0.00092(58) 


0.76 


1 




50 




0.972 




100 


4100 




0.13232 


-0.00079(38) 


-0.00109(38^ 


0.77 


1 




50 




0.947 


% 


100 


4000 
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/3 = 2.60, L/a = 8 



1.20 





13531 


0.02110 


64; 




0.00168(38; 




0.878 




"64 


0.979 


3) 


110 


4500 







13550 


0.01528 


44^ 




0.00142(85; 




0.879 




64 


0.972 


3 


110 


4500 







13574 


0.00810 


69; 




0.00170(43; 




0.87(1) [641 


0.966 


6) 


120 


4500 







13594 


0.00140 


72 




0.00158(62 




0.882 




64 


0.965 


8) 


130 


4500 


1.27 





13454 


0.02061 


53; 




0.00073(91, 




0.870 


?i 


64 


0.983 


2) 


110 


4500 







13473 


0.01512 


73; 




0.00218(47; 




0.881(6) 


64 


0.978 


2) 


110 


4500 







13496 


0.00721 


52^ 




0.00102(48; 




0.883(10) [641 0.971 


5 


120 


4500 







13516 


0.00327 


73; 




0.00039(53' 




0.883 


'6) 


64 


0.972 


3) 


130 


4500 


1.34 





13378 


0.02177 


75 




-0.00006(57; 


0.883 


6) 


64 


0.984 


2) 


110 


4500 







13420 


0.00830 


57^ 




-0.00004(45; 


0.872 


6) 


64 


0.972 


3) 


120 


4500 







13440 


0.0018(11) 




0.00031(44) 


0.876 


'9) 


64 


0.975 


2) 


130 


4500 







13473 


-0.00968(93) 


-0.00007(34; 


0.880 


6) 


64 


0.903 


5) 


110 


4500 


1.41 





13303 


0.02040(50' 




-0.00055(43; 


0.874 


7) 


64 


0.983 


2) 


110 


4500 







13322 


0.01375 


64; 




-0.00042(67; 


0.874 


'7) 


64 


0.979 


3) 


110 


4500 







13344 


0.00792 


81 




-0.00064(37 


0.873 


9) 


64 


0.975 


3) 


120 


4500 







13364 


0.00116 


66^ 




-0.00022(50, 


0.882 


6) 


64 


0.964 


3) 


130 


4500 


1.48 





13277 


0.0036(10) 




-0.00109(49; 


0.883 


'6) 


64 


0.978 


3) 


130 


5000 







13301 


-0.00316(79) 


-0.00085(34 


0.872 


6) 


64 


0.979 


2) 


150 


5000 


1.55 





13202 


0.00476(58) 


-0.00218(36; 


0.873 


7) 


64 


0.980 


2) 


130 


5000 







13226 


-0.00270(73) 


-0.00195(48' 


0.879(10) [64] 0.973(3) 


140 


5000 



1.3 



1.4 



1.5 



1.6 
1.7 



0.135917 
0.136152 
0.136387 
0.134882 
0.135113 
0.133410 
0.133636 
0.133862 
0.132400 
0.132680 
0.131230 
0.131510 



0.0211(39) 
0.01146(50' 
0.00276(49 
0.01207(50 
0.00466(45 
0.0206(10) 
0.01300(56 
0.00584(49 
0.0151(12) 
0.00567(70) 
0.01385(71) 
0.0047(13) 



P = 2.40, L/a = 8 

0.00162(74) 
0.00011(100) 

0.00100(29) 

0.00022(60 

0.00034(50 
-0.00040(30) 

0.00027(68) 
-0.00095(39) 
-0.00210(81' 
-0.00194(55 
-0.00236(54 
-0.00321(28 



0.819 




50 


0.974 


2) 


110 


10000 


0.819 


?i 


50 


0.971 


2) 


120 


10000 


0.82(1) 


50] 


0.947 


2) 


120 


10000 


0.815 


5) 


50 


0.974 


2) 


120 


10000 


0.828 


5) 


■50 


0.954 


2) 


120 


10000 


0.823 


4) 


50 


0.983 


2) 


110 


10000 


0.828 




'50' 


0.979 




120 


10000 


0.827 


•l\ 


'50' 


0.966 


I] 


120 


10000 


0.82(1) [501 


0.980 




120 


11900 


0.825(5) 


50 


0.966 


li 


120 


11900 


0.888 




64 


0.983 




120 


10700 


0.886 


'11 


64 


0.978 


5i 


130 


10700 



(3 = 2.20, L/a = 8 



1.3 


0.138247 


0.01685(90) 0.00161(63) 


0.840(5) 


50 


0.958(2) 


120 


16500 




0.138487 


0.0100(15) 


0.00144(36) 


0.836 




50 


0.946 




130 


16500 




0.138729 


0.00128(62' 


0.00149(45) 


0.843 


■li 


50 


0.919 




140 


16500 


1.5 


0.135400 


0.01877 


64; 


-0.00028(38; 


0.844 


'4) 


50 


0.973 


1) 


120 


16000 




0.135654 


0.01083 


46 


-0.00004(36; 


0.844 


4) 


50 


0.965 


2) 


130 


16500 




0.135885 


0.00285 


55, 


0.00079(51) 


0.841 


4) 


50 


0.951 


2) 


140 


16500 


1.7 


0.132712 


0.01913 


70; 


-0.00148(29' 


0.846 


'5) 


50 


0.983 


1) 


120 


16500 




0.132934 


0.01226 


45 


-0.00198(45 


0.844 


8) 


50 


0.977 


2) 


130 


16500 




0.133156 


0.00494 


49, 


-0.00205(50, 


0.834 


.5) 


50 


0.968 


2) 


140 


16500 


1.9 


0.130170 


0.02045 


54; 


-0.00385(32; 


0.843 




50 


0.984 


1) 


120 


16100 




0.130370 


0.01248 


76 


-0.00342(30 


0.840 


1 


50 


0.984 


2) 


130 


16100 




0.130570 


0.00561 


55 


-0.00322(32^ 


0.841 




50 


0.977 


1) 


140 


16100 



/3 = 2.10, L/a = 8 



1.5 





1355 


0.0579(12) 


0.00066(57) 


0.850(4) 


50 







1358 


0.0500(15) 


0.00011 


44) 


0.847 


4) 


50 







1360 


0.04380(77) 0.00033 


44) 


0.855 


4) 


■50 







1362 


0.0386(10) 


0.00047 


56) 


0.851 


.6) 


■50 


1.6 





1340 


0.06167 


75; 


-0.00027(49) 0.854 


'5) 


50 







1344 


0.04811 


59 


0.00039(83) 


0.848 


5) 


■50 







1346 


0.04179 


72 


-0.00035 




0.850 


4) 


■50 


1.7 





1326 


0.05964 


58, 


-0.00132 




0.853 


.5) 


■50 







1329 


0.04952 


50; 


-0.00128 


'33; 


0.850 


'4) 


'50 







1331 


0.04393 


48 


-0.00145 


32 


0.854 


3) 


■50 







1333 


0.03752 


54, 


-0.00018 


.31; 


0.849 


.3) 


■50 







1335 


0.02976 




-0.00096 




0.856 




'50 


1.8 





1315 


0.05054 


s: 


-0.00192 


ISi 


0.850 


3! 


■50 







1318 


0.04101 




-0.00182 


•s- 


0.844 




'50 







1321 


0.03149 




-0.00176 




0.847 


•l\ 


■50 







1324 


0.01956(59' 


-0.00148(33' 


0.850(4) 


■50 



0.979(3) [80] 
0.9877(10) [100] 

0.984(1) 100] 

0.978(2) [100] 

0.986(2) [80] 
0.9888(10) [100] 

0.985(1) [100] 

0.987(1)'""' 

0.987(2) 

0.988(1' 

0.982(2 

0.974(1 

0.9889(9) 
0.9883(10) 

0.983(2) 

0.980(1) 



90 
100' 
100' 
100' 
90' 
100 
110 
120' 



14500 
14500 
14400 
14500 
14500 
14500 
14500 
14500 
14500 
14500 
14500 
14500 
14500 
14500 
14500 
14500 



13 = 2.00, L/a = 8 



1.5 0.1383550 
0.1386672 



0.0276(15) 
0.0176(12) 



0.00020(44) 
0.00082(56) 



0.863(3) 
0.857(3) 



0.961(2) 
0.928(2) 



130 
140 



24500 
21500 



15 








1388000 


0.0138( 


13) 


0.00034(58) 





853 




50 


0.914(7)[160 


21200 


1.6 





1364310 


0.0360( 


10 


0.00036(60) 





858 


;4) 


50 


0.9912 


(9)fl40l 20000 







1367500 


0.02618(98) -0.0011(10) 





857 


(4) 


50 


0.982 




150 


20000 







1370700 


0.0156(12) 


0.00049(52) 





857 


7) 


50 


0.961 




160 


20000 


1.7 





1348740 


0.03723(61 


-0.00134(4i; 





856 


.3) 


50 


0.9932 


(9) [1401 20000 







1354309 


0.01588 


70 


-0.00073(30J 





861 


(3) 


50 


0.950 


2) 


130 


24500 







1354980 


0.01435 


93 


-0.00037(46 





859 


3) 


50 


0.971 


1) 


160 


20000 







1356000 


0.01043 


(82 


-0.00072(58; 





861 


4) 


50 


0.951 


2) 


150 


20300 


1.8 





1333520 


0.03533 


61 


-0.00095(40; 





859 


(3) 


50 


0.9936 




[140 


1 20000 







1336570 


0.02441 


57 


-0.00129(33 





863 


3) 


50 


0.9899 


if 


[150 


] 20000 







1339620 


0.01374 


(65 


-0.00050(51, 





853 


5) 


50 


0.978(1) 


160 


20000 


1.9 





1320578 


0.02843 


57 


-0.00237(38; 





860 


'3) 


50 


0.977 


1) 


110 


24500 







1323422 


0.01918 


68 


-0.00166(42^ 





860 


2) 


50 


0.976 


1) 


130 


24500 







1326278 


0.01016 


60 


-0.00183(39, 





858 


3) 


50 


0.945 


2) 


130 


24500 


2.0 





1308300 


0.0208(13) 


-0.00357(75; 





853 




50 


0.989 




150 


5400 







1311100 


0.01208(82 


-0.00310(40; 





860 


?i 


50 


0.983 




160 


7300 


2.1 





1293800 


0.02589 


(84 


-0.00391(65; 





859 




50 


0.996 




150 


5400 







1296600 


0.01483 


[96 


-0.00406(79^ 





855 


?i 


50 


0.988 




160 


7300 



P = 2.00, L/a = 6 



1 ^0 

-L .OU 


n 


1 Ann 


0.1002 


43; 


0.0055(15) 







1 


OyJ 


0.9998 


(2) 


120 


11200 




n 


1 4nn 


0.1017 


33 


0.00346( 


89) 









OyJ 


1.0000 


0) 


120 


11200 







1405 


0.0871 


30; 


0.00390( 


95j 





923 


2) 


50 


0.9993 




120 


14500 







1405 


0.0890 


29; 


0.0030( 







926 




50 


0.9991 




120 


14500 







1410 


0.0669 


28< 


0.0033( 


!!! 





922 


?i 


50 


0.9974 




120 


15000 







1410 


0.0721 


30; 


0.0045( 


15) 





922 


'2) 


50 


0.9974(8) 


120 


15000 







1415 


0.0505 


26 


0.0037( 


11) 





922 


*4) 


50 


0.990( 




120] 


15100 







1415 


0.0508 


24; 


0.0048( 


12) 





923 


2) 


50 


0.990( 


!! 


120] 


15100 


1.45 





1380 


0.0752 


22; 


0.00153(95) 





925 


'3) 


50 


0.9995 


(2) 


120 


14900 







1380 


0.0784 


17 


0.00126(66) 





925 


3) 


50 


0.9994 


3) 


120 


14900 







1385 


0.0564 


17; 


0.00224(73) 





921 


.3) 


50 


0.9985 


4) 


120 


14900 







1385 


0.0568 


24; 


0.00094(76) 





922 


'4) 


50 


0.9979(5) 


120 


14900 







1390 


0.0401 


20 


-0.0001(12; 





921 


2) 


50 


0.994( 




120 


15000 







1390 


0.0417 


21; 


0.00173(83) 





921 


3) 


50 


0.993( 


!! 


120 


15000 







1395 


0.0231 




0.00128(80) 





923 




50 


0.982( 




120 


15100 







1395 


0.0245 


23] 


0.00238(83) 





926 




50 


0.981( 


?! 


120 


15100 


1.60 





1355 


0.0673 


20; 


0.0021(16) 





925 


'4) 


50 


0.9994 


(2) 


120 


11200 







1355 


0.0698 


16 


-0.0009C 


(78 


} ° 


924 


3) 


50 


0.9998 


2) 


120 


11200 







1360 


0.0530 


17; 


-0.00053 


(56 





920 


3) 


50 


0.9983 


(6) 


120 


15000 







1360 


0.0535 


14; 


-0.0004(14) 


925 


'3) 


50 


0.9990 


3) 


120 


15000 







1365 


0.0340 


19 


0.0008(10) 





926 


2) 


50 


0.9957 


8) 


120 


15000 







1365 


0.0347 


15; 


0.00078(62) 





925 


3) 


50 


0.9971 


(5) 


120 


15000 







1370 


0.0165 


14; 


-0.00013(69) 


922 


'3) 


50 


0.988( 


1) 


120] 


15000 







1370 


0.0171 


16 


-0.0004(12) 


923 


2) 


50 


0.9888( 


10) [1201 15000 


1.75 





1330 


0.0682 


37 


-0.00223(95) 


920 


7) 


50 


1.0000(0) 


120 


3400 







1330 


0.0683(20) 


-0.0016(10) 


927 


7) 


50 


0.9997 


(3) 


120 


3400 







1335 


0.05035 


(96 


) -0.00293 


(66 


) 


925 


3) 


50 


0.9996 


3) 


120 


14500 







1335 


0.05242 


(96 


) -0.00187 


(64 


) 


923 


3) 


50 


0.9995 


2) 


120 


14500 







1340 


0.0339(10' 


-0.0008S 


(59 


) 


927 


.3) 


50 


0.9985 


(5) 


120 


15000 







1340 


0.0350 




-0.0025 


[!? 





924 




50 


0.9983 




120 


15000 







1345 


0.0164 




-0.0026 







920 


:3i 


50 


0.9948 




120 


15000 







1345 


0.0197(15] 


-0.00095(70 


) 


920 


^3) 


50 


0.994(1) [1201 


15000 



/3 = 1.90, L/a = 6 






1410 


0.1315( 


621 


0.00112(92) 





930 


'2) 


50 


0.9997 


'1) 


120 


24000 





1410 


0.1338( 


41 


0.0008(13) 





923 


2) 


50 


0.9997 


1) 


120 


23900 





1415 


0.1103( 


64 


0.0006(12) 





928 


2) 


50 


0.9985 


4) 


120 


24000 





1415 


0.1136( 


36 


0.00144(94) 





927 




50 


0.9989 




120 


24000 





1420 


0.0857 


25) 


0.0028(10) 





930 


?i 


50 


0.9930 


?! 


120 


24100 





1420 


0.0876 


31; 


0.00166(90) 





923 


'2) 


50 


0.9945 


'5) 


120 


24100 





1340 


0.0878( 


25 


-0.00021 







929 


3) 


50 


1.0000 


0) 


120 


15100 





1340 


0.0900 


21; 


-0.00174 


lis! 





928 


.3) 


50 


1.0000 


.0) 


120 


15100 





1345 


0.0710 


21; 


-0.00159 


'69; 





926 


'3) 


50 


0.9994 


'2) 


120 


24000 





1345 


0.0733 


14 


-0.00193 


55 





929 


2) 


50 


0.9995 


3) 


120 


24000 





1350 


0.0520 


24; 


-0.00136 


.59; 





929 


2) 


50 


0.9973 


.3) 


120 


24100 





1350 


0.0535 




-0.00045 







928 




50 


0.9976 




120 


24100 





1355 


0.0307 


Si 


-0.00102 







926 


?! 


■50 


0.9888 




120 


24100 





1355 


0.0309(14) 


-0.00108 







926 




'50 


0.9895(8) 


120 


24100 





1280 


0.06146(84 


) -0.00601 







928 


?i 


50 


0.99991(9) [120] 22700 





1280 


0.06303 


(92 


) -0.00699 







926 




50 


0.9998 




120 


22700 





1285 


0.04376 


(79 


) -0.00734 


•Si 





928 


■ii 


50 


0.9993 


?i 


120 


23600 





1285 


0.04434 


(83 


) -0.00594 







928 




50 


0.9996 




120 


23600 





1290 


0.02540 


(87 


) -0.00681 







929 


§i 


50 


0.9976 




120 


24100 



0.1290 


0.02607(99) 


-0.00585(67) 0.924 


>l\ 


50 


0.9976(4) 


120 


24100 


n 1 one; 


U.UU < 4/ (99 ) 


-0.0044(16) 


0.925 


^) 


OU 


0.9916(8) 


120 


24100 


0.1295 


0.0084(15) 


-0.00649(46^ 


0.928 


'2) 


50 


0.9922(6) 


120 


24100 


2.5 0.1240 


0.0535(12) 


-0.01090(60 


0.926 


2) 


50 


0.99990(7) [1201 22200 


0.1240 


0.0539(12) 


-0.01076(50; 


0.927 


2) 


50 


0.9997(2) 


120 


22200 


0.1245 


0.0372(13) 


-0.01083(49^ 


0.929 


'2) 


50 


0.9996(1) 


120 


24000 


0.1250 


0.01913(94) 


-0.01152(87 


0.926 


2) 


50 


0.9985(3) 


120 


24000 


0.1250 


0.0198(11) 


-0.01144(38* 


0.925 


2) 


50 


0.9986(3) 


120 


24100 


0.1255 


0.00201(98) 


-0.01160(49) 


0.924 




50 


0.9937(7) 


120 


24100 


0.1255 


-0.00014(87) 


-0.01100(44; 


0.923 


■1 


50 


0.9932(8) 


120 


24100 



TABLE VI: Same as Table|V|for two-flavor QCD. 



Csw 




aM 


aAM 


acc 


|iVMD 


iVtraj 






(3 = 12.00, L/a = 


8 










1.00 


0.12659 


0.01266(15) 


0.00082(16) 


0.75(3) 


100 


1500 




0.12676 


0.007137(77) 


0.00083( 


14^ 


0.74 




100 


1500 




0.12693 


0.00180(16) 


0.00115( 


16) 


0.73 


■ii 


100 


1500 




0.12709 


-0.003231(98) 0.00070( 


13) 


0.75 


2) 


100 


1500 


1.05 


0.12659 


0.00882(14) 


0.00003( 


20) 


0.71 


3) 


100 


1500 




0.12676 


0.00317(11) 


0.00019( 


14) 


0.73 


2) 


100 


1500 




0.12693 


-0.00199(11) 


0.00016(14) 


0.70 


'2) 


100 


1500 




0.12709 


-0.00735(16; 


-0.00008 


(12 


) 0.72 


2) 


100 


1500 


1.10 


0.12659 


0.00498(12) 


-0.00061 


(11 


) 0.74 


2) 


100 


1500 




0.12676 


-0.00042(12' 


-0.00065 


(15 


) 0.73 


1) 


100 


1500 




0.12693 


-0.00581(11 


-0.00069 


(12 


) 0.71 


2) 


100 


1500 




0.12709 


-0.01119(15; 


-0.00070 


(12 


) 0.74 


2) 


100 


1500 



/3 = 8.85, L/a = 8 



1.040 


0.1270 


0.020524(91) 


0.00036(11) 


0.72(2) 


80 


2100 




0.1274 


0.00794(13) 


0.00056 


;ii; 


0.72 


(2) 


80 


2100 




0.1278 


-0.00470(14) 


0.00052 


15 


0.70 


1) 


80 


2100 




0.1282 


-0.017306(94) 


0.00043 


;i4; 


0.70 


;i) 


80 


2000 


1.055 


0.1270 


0.019165(83) 


0.00044 


'14) 


0.70 


'2) 


80 


2000 




0.1274 


0.00665(14) 


0.00042 


16 


0.70 


2) 


80 


2000 




0.1278 


-0.00628(12) 


0.00032 


;i8; 


0.67 


;2) 


80 


2000 




0.1282 


-0.01900(12) 


0.00025 


'16) 


0.70 


(2) 


80 


2000 


1.070 


0.1270 


0.01758(10) 


0.00002 


15 


0.70 


1) 


80 


2000 




0.1274 


0.00505(12) 


0.00001 


18 


0.70 


1) 


80 


2000 




0.1278 


-0.00754(14) 


0.00027 


'14* 


0.69 


2) 


80 


2000 




0.1282 


-0.02057(17) 


0.00012(11' 


0.70(2) 


80 


2000 



/3 = 5.00, Lja = 8 



1.09 


0.12954 


0.01204(16) 


0.00054( 


26) 


0.75 


;i) 


64 


2300 




0.12970 


0.00692(23) 


0.00068( 


21) 


0.76 


2) 


64' 


2300 




0.12986 


0.00198(25) 


0.0004l( 


21) 


0.74 


.1) 


64' 


2300 




0.13002 


-0.00308(18) 


0.00086(43) 


0.74 


;i) 


64' 


2300 


1.13 


0.12933 


0.01167(21) 


-0.00005 


(31 


) 0.74 


2) 


64 


2300 




0.12949 


0.00691(21) 


-0.00014 


(24 


) 0.75 


.1) 


64 


2300 




0.12965 


0.00174(15) 


0.00041(19) 


0.75 


;i) 


64 


2300 




0.12981 


-0.00307(24) 


-0.00013 


(20 


) 0.75 


1) 


64 


2300 


1.17 


0.12912 


0.01164(26) 


-0.00063 


(18 


) 0.76 


.2) 


64 


2300 




0.12928 


0.00667(27) 


-0.00031 


(28 


) 0.74 


'2) 


64' 


2300 




0.12943 


0.00175(32) 


-0.00076 


(20 


) 0.75 


1) 


64 


2300 




0.12959 


-0.00338(19) 


-0.00020 


(33 


) 0.75 


1) 


64 


2300 



/3 = 3.00, Lja = 8 



1.20 


0.1332100 


0.02195 




0.00069(38) 


0.780(5) 


50 


10500 




0.1333700 


0.01674 




0.00086(22) 


0.767(7) 


50 


10500 




0.1335400 


0.01165 


'34; 


0.00044(28) 


0.77(1) [501 


10500 




0.1337000 


0.00685 


'22* 


0.00068(28) 


0.762(5) 


50 


10500 


1.28 


0.1324700 


0.02175 


35; 


-0.00013(29 


) 0.774(6) 


50 


10500 




0.1326300 


0.01668 




-0.00007(21 


) 0.773(6) 


50 


10500 




0.1327900 


0.01168 




-0.00008(40 


) 0.781(8) 


50 


10500 




0.1329500 


0.00649 


'9 


0.00021(22) 


0.766(5) 


50 


10500 


1.36 


0.1317300 


0.02146 




-0.00079(21 


) 0.769(7) 


■50 


10500 




0.1318900 


0.01644 


'30; 


-0.00068(24 


) 0.771(7) 


50 


10500 




0.1320400 


0.01156 


37 


-0.00099(28 


) 0.777(7) 


50 


10500 




0.1322000 


0.00675 


30* 


-0.00092(22 


) 0.779(6) 


■50 


10500 
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TABLE VI: {Continued) 



P = 2.60, L/a = 8 



1.20 



1.25 



1.30 



1.35 



0.135574 
0.135738 
0.135903 
0.136068 
0.135020 
0.135180 
0.135340 
0.135510 
0.134470 
0.134630 
0.134790 
0.134950 
0.133920 
0.134080 
0.134240 
0.134400 



0.02473(74) 
0.0191(12) 
0.01522(75) 
0.01002(58' 
0.02424(81 
0.01912(58 
0.01506(67 
0.00982(59) 
0.02492(54) 
0.01950(73) 
0.01474(63) 
0.00949(53) 
0.02369(45) 
0.0190(10) 
0.01426(70) 
0.00907(51) 



0.00104(62) 
0.00161(59) 
0.00209(53) 
0.00143(93) 
0.0015(10) 
0.00149(49' 
0.00156(56 
0.00150(82) 
-0.0004(11) 
0.00175(52) 

-0.00007(49) 
0.00038(65) 
0.00044(47) 

-0.00040(73) 
0.00011(93) 

-0.00044(57) 



0.809(7) 
0.800(8) 
0.808(7) 
0.806(8) 



0.811(8 
0.810(8' 
0.812(7 
0.800(8 
0.814(8' 
0.803(7 
0.804(8' 
0.82(1)[50' 
0.803(7) [50] 
0.82(1)' ■ 
0.81(1) 



4500 
4500 
4500 
4500 
4500 
4500 
4500 
4500 
4500 
4500 
4500 
4500 
4500 
4500 
4500 
4500 



(3 = 2.20, L/a = : 



1.35 


0.13868 


0.01252(59) 


0.00182(70' 


0.828(3) 


50 


37300 




0.13914 


-0.00310(93) 


0.00130 


87; 


0.834 


2) 


50 


35900 


1.50 


0.13654 


0.00838(96) 


0.00071 


72 


0.834 


*2) 


50 


41100 




0.13693 


-0.00433(75) 


0.00003 


62^ 


0.833 


2) 


50 


39600 


1.60 


0.13500 


0.01012(23) 


-0.0006C 


)(39) 0.839 




50 


39900 




0.13543 


-0.00483(38) 


-0.00076(68) 0.832 




50 


38000 



1.38 



1.53 



1.63 



1.73 



0.14040 
0.14092 
0.13741 
0.13837 
0.13599 
0.13648 
0.13451 
0.13497 



(3 = 2.10, L/a = 8 
0.00598(95) 



-0.01094(54) 
0.02124(23) 

-0.00872(53) 
0.01371(17) 
-0.0018(11) 
0.01152(32) 

-0.00317(40) 



0.0007(11) 
0.00188(56) 
0.00062(65) 
0.00074(92) 
0.00030(58) 
0.00020(57) 
-0.00008(48) 
-0.0018(12) 



0.811(2) 
0.810(2) 
0.829(2) 
0.822(2) 
0.830 
0.828 
0.836(1 
0.833(2 



145700 
167500 
67800 
58900 
66400 
62900 
104400 
137100 



1.2 



1.4 



1.6 



1.8 



2.0 



2.4 



0.14347 
0.14391 
0.13987 
0.14021 
0.14056 
0.13660 
0.13683 
0.13725 
0.13335 
0.13362 
0.13389 
0.13059 
0.13090 
0.12500 
0.12550 



f3 = 2.10 

0.0077(22) 
-0.0025(21) 

0.0125(18) 

0.0043(14) 
-0.0090(23) 

0.0073(13) 
0.00065(93) 
-0.0136(16) 
0.00825(67) 
-0.00169(79) 
-0.0102(10) 
-0.0001(14 
-0.0127(12 

0.0022(17) 
-0.0148(31) 



, L/a = 6 

0.0063(21) 
0.0079(24) 
0.0030(11) 
0.0032(17) 
0.0025(14) 
0.0014(10) 
-0.00042(78) 
-0.0002(11) 
-0.00438(74) 
-0.00214(84) 
-0.00278(77) 
-0.00685(74) 
-0.0071(10) 
-0.01283(75) 
-0.0133(12) 



0.864(3) 
0.865(4) 
0.868(4) 
0.864(3) 
0.862(4) 
0.869(7) 
0.868(3) 
0.863(2) 
0.867(4) 
0.870(2) 
0.867(4) 
0.870(3) 
0.870(3) 
0.873(3) 
0.865(3) 



21600 
26000 
26000 
26000 
26000 
25200 
26000 
25200 
26000 
26000 
28000 
21600 
21600 
21600 
21600 



1.4 



1.6 



1.8 



2.0 



2.2 



2.6 



13 = 2.00, L/a = 6 



0.14279 
0.14362 
0.13901 
0.13936 
0.13971 
0.13536 
0.13588 
0.13627 
0.13192 
0.13221 
0.13250 
0.12839 
0.12902 
0.12295 
0.12353 



0.0125(26) 
-0.0135(21) 

0.0131(17) 

0.0035(16) 
-0.0088(16) 

0.0113(15) 
-0.0069(12) 
-0.0198(45) 

0.0116(11) 

0.0015(18) 
-0.0106(12) 

0.0207(11) 
-0.0008(18) 
0.01399(86) 
-0.0093(14) 



0.0068(17) 
0.0033(33) 
0.0034(21) 
0.00330(91) 
0.0009(26) 
0.0002(13) 
-0.0013(22) 
0.0006(24) 
-0.00490(92) 
-0.0041(14' 
-0.0033(11 
-0.0074(11 
-0.0078(12 
-0.01330(73 
-0.01204(69 



0.906( 
0.902 ( 
0.914( 
0.907( 
0.906( 
0.918( 
0.909( 
0.907( 
0.918( 
0.916( 
0.915( 
0.921( 
0.919( 
0.918( 
0.921( 



20000 
25000 
25000 
25000 
25000 
25000 
25000 
25000 
25000 
25000 
25000 
20000 
20000 
21000 
21000 



18 



TABLE VI: (Continued) 



TABLE VII: Same as Table |V| for quenched QCD. 



csw 




aM 


aAM 




-face 


|A'"md| 


iVtraj 






P = 24.00, L/a = 


= 8 










1.00 


0.12567 


0.009136(59) 


0.000714( 


97) 


0.67 




128 


3100 




0.12584 


0.003733(51) 


0.000545( 


69) 


0.67 


l!i 


128 


3100 




0.12600 


-0.00139(11) 


0.000525( 


96) 


0.65 


'2) 


128 


3100 




0.12617 


-0.006783(76) 


0.000567( 


79) 


0.66 


r\ 


128 


3100 


1.03 


0.12567 


0.007822(63) 


0. 000078 ( 


90) 


0.67 




128 


3100 




0.12584 


0.002319(54) 


0.000054( 


70) 


0.67 


1 


128 


3100 




0.12600 


-0.002756(68) 


0.000167( 


81) 


0.67 




128 


3100 




0.12617 


-0.008282(63) 


-0.000028 


(77 


) 0.67 




128 


3100 


1.06 


0.12567 


0.006427(46) 


-0.000334 


(63 


) 0.67 


ili 


128 


3100 




0.12584 


0.000984(78) 


-0.000317 


(76 


) 0.68 


;i) 


128 


3100 




0.12600 


-0.00401(10) 


-0.000293 


(75 


) 0.65 


2) 


128 


3100 




0.12617 


-0.009544(91) 


-0.000295 


(78 


) 0.66 


1) 


128 


3100 



P = 12.00, L/a = 8 



1.00 


0.12659 


0.013236(96) 


0.00080( 




0.74 




100 


1600 




0.12676 


0.007870(89) 


0.00086( 




0.73 




100 


1600 




0.12693 


0.00246(14) 


0.00086( 


17) 


0.72 


'1) 


100 


1600 




0.12709 


-0.00261(10) 


0.00107( 


15) 


0.73 


2) 


100 


1600 


1.05 


0.12659 


0.009444(92) 


0.00023( 


14) 


0.74 


2) 


100 


1600 




0.12676 


0.004108(91) 


0.00012(12) 


0.72 


2) 


100 


1600 




0.12693 


-0.001424(88) 


-0.00000(15) 


0.70 


2) 


100 


1600 




0.12709 


-0.006318(91) 


0.000214(98) 


0.72 


[2) 


100 


1600 


1.10 


0.12659 


0.005635(81) 


-0.00043 


(13) 


0.73 


[2) 


100 


1600 




0.12676 


0.00026(16) 


-0.00048 


(11) 


0.74 


1) 


100 


1600 




0.12693 


-0.00507(15) 


-0.00053 


(10) 


0.73 


2) 


100 


1600 




0.12709 


-0.01021(11) 


-0.00069(16) 


0.75(2) 


100 


1600 



(3 = 8.85, L/a = 8 



1.05 


0.1261 


0.048752(68) 


0.000670(83) 


0.70 


1) 


80 


3500 




0.1266 


0.033381(77 


0.000472(91 


0.70 


1) 


80 


3500 




0.1271 


0.017522(76) 


0.000515(90) 


0.69 


1) 


80 


3500 




0.1276 


0.001942(99' 


0.00034(11) 


0.69 


1) 


80 


3500 


1.07 


0.1260 


0.049789(63 


0.000197(85) 


0.70 


l) 


80 


3500 




0.1265 


0.034276(80 


0.000152(88) 


0.68 


1) 


80 


3500 




0.1270 


0.018797(76) 


0.00027(12) 


0.70 


1) 


80 


3500 




0.1275 


0.00324(11) 


0.00019(12) 


0.69 


1) 


80 


3500 


1.09 


0.1259 


0.051029(60) 


-0.000008(97) 0.69 




80 


3500 




0.1264 


0.035509(77; 


-0.000033(81) 0.70 




80 


3500 




0.1269 


0.019930(91) 


-0.00006(12) 0.71 




80 


3500 




0.1274 


0.00427(13) 


0.00005(12) 


0.70 


■!i 


80 


3500 









/3 = 





12954 


0.01759 







12970 


0.01254 







12979 


0.00991 


18) 





12986 


0.00767 


17) 





12995 


0.00480 


21) 





13002 


0.00280(14) 





13011 


-0.00014 


(16) 





13027 


-0.00508(20) 



1.13 0.12951 
0.12954 
0.12967 



0.01003(15) 
0.009028(95) 






12970 


0.00413( 


13) 





12983 


0.00026( 


12) 





12986 


-0.00082 


(15 





12999 


-0.00489 


(18 





13002 


-0.00576 


(11 





12924 


0.00995( 


12) 





12940 


0.00500( 


16) 





12954 


0.00075( 


15) 



0.12956 
0.12970 



0.00010(14) 
-0.00429(23) 



(3 = 5.00, L/a = 8 

0.00080(14' 
0.00075(25 
0.00057(17) 
0.00100(16) 
0.00115(19) 
0.00083(25) 
0.00094(27) 
0.00069(27) 
-0.00002(18) 
0.00041(12' 
0.00012(26* 
-0.00009(29) 
0.00048(18) 
0.00006(25) 
0.00038(23) 
0.00019(20) 
-0.00052(19) 
-0.00023(25) 
-0.00012(21) 
-0.00005(20) 
-0.00019(32) 



0.762(10) [64] 
0.751(9) [64] 
0.73(1) [64] 
0.754(9) [64] 
0.75(1)' ' 
0.75(1) 



64 
64 

0.738(9)'[64] 
■ ' '64' 
64 



0.76(1) 
0.75(2) 
0.747(9)'[64 
0.75(1) [64]" 
0.725(10) [64] 
0.74(1) [64] 
0.75(2) 64 
0.733(10) [64] 
0.747(9) [64] 
0.74(1) [64] 
0.75(2) 64] 
0.742(10) [64] 
0.75(1) [64' 
0.73(2) 64 



4800 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
3500 
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TABLE VII: {Continued) 



0.12972 


-0.00515 


16) 


-0.00054 


22) 


0.75(1) 


■64] 


3500 


0.12986 


-0.00933 


13) 


-0.00018 


27) 


0.74(2) 


64 


3500 


0.13002 


-0.01466 


12) 


-0.00032 


17) 


0.755(9) [64] 


3500 



(3 = 3.00, L/a = 8 



1.20 





13393 


0.01524(24) 


0.00083(26) 




0.785(6) 


50 


8100 







13410 


0.01010 


25 


0.00144(32) 




0.789 


6) 


50 


8100 







13428 


0.00461 


25 


0.00152(41 




0.786 


8) 


50 


8100 







13440 


0.00120 


27) 


0.00132(59) 




0.778 


6) 


50 


8100 


1.28 





13315 


0.01562 


28) 


0.00057(29) 




0.77(1) [50] 


8100 







13332 


0.01053 


18 


0.00010(29 




0.782 




50 


8100 







13349 


0.00520 


3o; 


0.00122(38) 




0.782 


?! 


50 


8100 







13363 


0.00110 


23) 


0.00003(36' 




0.777 




50 


8100 


1.36 





13239 


0.01422 


36; 


-0.00045(27) 


0.777 




50 


8100 







13255 


0.00941 




-0.00050(34) 


0.776 


7) 


50 


8100 







13272 


0.00414( 




0.00008(33) 


0.777 


8) 


50 


8100 







13286 


-0.00033(36) -0.00049(35) 


0.776 


9) 


50 


8100 



/? = 2.70, L/a = 8 



1.2 





13605 


0.01302(55) 


0.00276( 


40) 







803 


4) 


50 


15000 







13642 


0.00235(76) 


0.00270( 


50^ 







803 


4) 


50 


15000 







13680 


-0.01013(71) 


0.0029(12) 







797 


4) 


50 


15000 


1.3 





13472 


0.01845(90) 


0.00187( 


60) 







798 


8) 


50 


15000 







13526 


0.00207(59) 


0.00158( 


62 







797 


4) 


50 


15000 







13544 


-0.00305(17) 


0.00150( 


24 







800 


5) 


50 


15000 


1.4 





13356 


0.01833(30) 


0.00015( 


50) 







805 


4) 


50 


15000 







13412 


0.00159(30) 


0.00113( 


48) 







798 


8) 


50 


15000 







13428 


-0.00405(34) 


0.00047( 


52< 







799 


7) 


50 


15000 


1.5 





13264 


0.01093(29) 


-0.00113(28) 





800 


5) 


50 


15000 







13281 


0.00562(27) 


-0.00107 


(37) 





805 


4) 


50 


15000 







13298 


-0.0009(12) 


-0.00076 


43 





799 


7) 


50 


15000 


1.6 





13143 


0.01254(23) 


-0.0037(14) 





797 


4) 


50 


15000 







13180 


0.00022(23) 


-0.00202(31) 





815(6) 


50 


15000 



(3 = 2.70, L/a = 6 



1.2 





13605 


0.01009(38) 


0.00323 




13 







888 


3) 


50 


15000 







13642 


-0.00071(39) 


0.00324 


48 







888 


3 


50 


15000 







13680 


-0.01069(42) 


0.00298 


64 







890 


4) 


50 


15000 


1.3 





13472 


0.01651(28) 


0.00106 


46 







890 


3) 


50 


15000 







13526 


0.00090(34) 


0.00097 


53 







890 


4) 


50 


15000 







13544 


-0.00499(36) 


0.00137 


44 







889 


6) 


50 


15000 


1.4 





13356 


0.01695(39) 


-0.00102 


) 


'36) 





889 


3 


50 


15000 







13412 


0.00017(32) 


-0.00088 


53 





890 


3) 


50 


15000 







13428 


-0.00464(39) 


-0.00070 


51 





890 


3 


50 


15000 


1.5 





13264 


0.00959(32) 


-0.00261 


41 





889 


6) 


50 


15000 







13298 


-0.00068(35) 


-0.00267 


*42 





890 




50 


15000 


1.6 





13143 


0.01083(32) 


-0.00504 


54 





890 


I] 


50 


15000 








,3 = 


2.60, L/a = 




6 















1.2 





13698 


0.00983(59) 


0.00469(88) 





896(3) 


50 


15000 







13730 


0.00064(65) 


0.00482(99) 





896 


3) 


50 


15000 







13749 


-0.00550(60) 


0.0049(11) 





893 


7) 


50 


15000 


1.3 





13574 


0.01057(49) 


0.00281(79) 





893 


7) 


50 


15000 







13616 


-0.00198(52) 


0.00170(81) 





894 




50 


15000 


1.4 





13463 


0.00750(57) 


-0.00051(53) 





894 


1 


50 


15000 







13494 


-0.00187(34) 


0.00068(58) 





894 


3) 


50 


15000 


1.5 





13331 


0.01048(54) 


-0.00249(50) 





891 


3) 


50 


15000 







13367 


-0.00005(38) 


-0.00176(56) 





893 


3) 


50 


15000 


1.6 





13215 


0.00774(45) 


-0.00342(50) 





894 




50 


15000 


1.8 





12953 


0.01242(62) 


-0.00778(50^ 





893 


1 


50 


15000 
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TABLE VIII: Numerical values of csw{go,L/a) and K{go,L/a) satisfying Eq. 13111 in three-flavor QCD. 



13 


L/a 


function 


csw(5'oi L/a) 


k(po, L/a) 


12.00 


8 


linear 


1.0546(25) 


0.1268421(61) 


8.85 


8 


linear 


1.0761(32) 


0.127513(10) 


5.00 


8 


linear 


1.1311(48) 


0.129641(26) 


3.00 


8 


linear 


1.254(15) 


0.13318(14) 


2.60 


8 


linear 


1.359(13) 


0.13423(14) 


2.40 


8 


linear 


1.384(23) 


0.13545(29) 


2.20 


8 


linear 


1.508(29) 


0.13587(39) 


2.10 


8 


linear 


1.649(58) 


0.13521(85) 


2.00 


8 


quadratic 


1.670(56) 


0.13639(89) 


2.00 


6 


quadratic 


1.632(45) 


0.13696(77) 


1.90 


6 


quadratic 


1.739(53) 


0.13741(98) 


TABLE IX: Same as Table IVTTTl but for two-flavor QCD. 


/3 


L/a 


function 


csw{go,L/a) 


K{g^,L/a) 


12.00 


8 


linear 


1.0558(27) 


0.1268509(66) 


8.85 


8 


linear 


1.0818(85) 


0.127519(27) 


5.00 


8 


linear 


1.1334(62) 


0.129686(34) 


3.00 


8 


linear 


1.276(20) 


0.13320(19) 


2.60 


8 


linear 


1.327(49) 


0.13496(55) 


2.20 


8 


linear 


1.519(32) 


0.13649(48) 


2.10 


8 


linear 


1.672(65) 


0.1358(10) 


2.10 


6 


linear 


1.598(19) 


0.13689(31) 


2.00 


6 


linear 


1.777(27) 


0.13612(47) 


TABLE X: Same as Table IVlTTI but for quenched QCD. 


f3 


L/a 


function 


csw{go,L/a) 


K{g^,L/a) 


24.00 


8 


linear 


1.0375(16) 


0.1259026(23) 


12.00 


8 


linear 


1.0627(27) 


0.1268574(65) 


8.85 


8 


linear 


1.0829(47) 


0.127565(15) 


5.00 


8 


linear 


1.1540(42) 


0.129701(23) 


3.00 


8 


linear 


1.338(20) 


0.13308(19) 


2.70 


8 


linear 


1.429(10) 


0.13380(12) 


2.70 


6 


linear 


1.3608(79) 


0.134554(92) 


2.60 


6 


linear 


1.414(14) 


0.13470(17) 
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T'A T>T T? VT. 

lAriLii- Ai: 


Final results for csw {go , L* / a) and K{go,L*/a) 


for fixed physical size L* for three-flavor QCD. 




L/a 


L*/a 


csw(q?), L* / a) 


Kc(an, L* /a) 


12.00 


8 


7.508095 xlO® 


1.0601(25) 


0.1269060(61) 


8.85 


8 


8.462365 xlO'* 


1.0903(32) 


0.127592(10) 


5.00 


8 


3.807760 xlO^ 


1.1825(48) 


0.129764(26) 


3.00 


8 


2.502040 xlO^ 


1.368(15) 


0.13325(14) 


2.60 


8 


1.475172 xlO^ 


1.467(13) 


0.13424(14) 


2.40 


8 


1.136512 xlO^ 


1.471(23) 


0.13544(29) 


2.20 


8 


8.780129 


1.542(29) 


0.13586(39) 


2.10 


8 


7.726477 


1.633(58) 


0.13521(85) 


2.00 


8 


6.805369 


1.576(56) 


0.13642(89) 


2.00 


6 


6.805369 


1.742(45) 


0.13691(77) 


2.00 


6.805369 


6.805369 


1.650(51) 


0.13669(83) 


1.90 


6 


6 


1.739(53) 


0.13741(98) 


TABLE XII: Same as Table IXll but for two-flavor QCD. 




L/a 


L*/a 


csw(ffo,-f'*/a) 


n{gl,L* /a) 


12.00 


8 


2.350129 xlO^ 


1.0583(27) 


0.1269100(66) 


8.85 


8 


3.656026 xW 


1.0919(85) 


0.127592(27) 


5.00 


8 


2.446546 xlO^ 


1.1797(62) 


0.129794(34) 


3.00 


8 


1.982120 xlO^ 


1.371(20) 


0.13323(19) 


2.60 


8 


1.219378 xlO^ 


1.405(49) 


0.13495(55) 


2.20 


8 


7.575548 


1.498(32) 


0.13651(48) 


2.10 


8 


6.738767 


1.590(65) 


0.1358(10) 


2.10 


6 


6.738767 


1.682(19) 


0.13680(31) 


2.10 


6.738767 


6.738767 


1.631(39) 


0.13638(63) 


2.00 


6 


6 


1.777(27) 


0.13612(47) 


TABLE XIII: Same as TableEH but for quenched QCD. 


P 


L/a 


L'/a 


csw(s'o,-f'*/a) 


n{gl,L*/a) 


24.00 


8 


3.088560 xlO" 


1.0264(16) 


0.1259370(23) 


12.00 


8 


2.409888 xlO'^ 


1.0590(27) 


0.1269080(65) 


8.85 


8 


6.326167 xlO^ 


1.0852(47) 


0.127627(15) 


5.00 


8 


8.042260 xlO^ 


1.1921(42) 


0.129763(23) 


3.00 


8 


9.115448 


1.359(20) 


0.13305(19) 


2.70 


8 


6.655769 


1.381(10) 


0.13389(12) 


2.70 


6 


6.655769 


1.4006(79) 


0.134470(92) 


2.70 


6.655769 


6.655769 


1.388(9) 


0.13426(10) 


2.60 


6 


6 


1.414(14) 


0.13470(17) 
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TABLE XIV: -L*/a, Scsw and Skc with the three- loop fJ function Eq. 1531 . 



p 


L j a 


L 1 Ob 


0CQ,v^{gQ,Lla] L j a) 


dnc\go, 1^/0', ^ /a) 


12.00 


8 


i.DooUD4 XiU 


5.5U9123X1U 


D.o489o2Xl(J 


8.85 


8 


1.755852 xlu 


1.412367x10 


7.951364x10 


5.00 


8 


o.44/o7o xlU 


4.79442UX lU 


I.z5l8l4x 10 


3.00 


8 


3.300322 xlO^ 


1.179115x10^^ 


9.547393x10"^ 


2.60 


8 


1.793750 xlO^ 


1.198918x10"^ 


3.523283x10"^ 


2.40 


8 


1.317591 xlO^ 


1.058685xl0~^ 


1.836828x10"'^ 


2.20 


8 


9.648020 


6.105866x10"^ 


-1.406025x10"^ 


2.10 


8 


8.244442 


1.257961x10"^ 


-4.125600x10"^ 


2.00 


8 


7.037491 


-7.046997x10"^ 


2.824179x10"^ 


2.00 


6 


7.037491 


1.335089x10"^ 


-6.098453x10"^ 


1.90 


6 


6 









23 



0.003 
0.002 
0.001 


-0.001 
-0.02 



p=12.0, L/a=8 



FIG. 1: aM dependence of oAM in three-flavor QCD. 
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FIG. 2: Same as Fig.[Tl but in two-flavor QCD. 
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FIG. 3: Same as Fig. but in quenched QCD. 
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